Multichannel quantum-defect theory for slow atomic collisions 
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We present a multichannel quantum-defect theory for slow atomic collisions that takes advan- 
tages of the analytic solutions for the long-range potential, and both the energy and the angular- 
momentum insensitivities of the short-range parameters. The theory provides an accurate and 
complete account of scattering processes, including shape and Feshbach resonances, in terms of a 
few parameters such as the singlet and the triplet scattering lengths. As an example, results for 
23 Na- 23 Na scattering are presented and compared close-coupling calculations. 
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I. INTRODUCTION 



Slow atomic collisions are at the very foundation of 
cold-atom physics, since they determine how atoms in- 
teract with each other and how this interaction might be 
manipulated yj, |2j ■ While substantial progress has been 
made over the past decade , there are still areas where 
the existing theoretical framework is less than optimal. 
For example, all existing numerical methods may have 
difficulty with numerical stability in treating ultracold 
collisions in partial waves other than the s wave, because 
the classically forbidden region grows infinitely wide as 
one approaches the threshold. This difficulty becomes 
a serious issue when there is a shape resonance right at 
or very close to the threshold, as the usual argument 
that the s wave scattering dominates would no longer be 
applicable. Another area where a more optimal formu- 
lation is desirable is analytic representation. Since much 
of our interest in cold atoms is in complex three-body 
and many-body physics, a simple, preferably analytical 
representation of cold collisions would not only be very 
helpful to experimentalists, but also make it much easier 
to incorporate accurate two-body physics in theories for 
three and many-atom systems. Existing formulations of 
cold collisions provide little analytical results especially 
in cases, such as the alkali-metal atoms, where the atomic 
interaction is complicated by hyperfine structures. Fur- 
thermore, whatever analytic results that we do have have 
been based almost exclusively on the effective-range the- 
ory the applicability of which is severely limited by 
the long-range atomic interaction 0, [|| . 

Built upon existing multichannel quantum-defect the- 
ories that are based either on free-particle reference func- 
tions or on numerical solutions for the long-range poten- 
tial 0, S H IH |H 113, we present here a multichan- 
nel, angular-momentum-insensitive, quantum-defect the- 
ory (MAQDT) that overcomes many of the limitations of 
existing formulations. It is a generalization of its single- 



channel counterpart and takes full advantage 

of both the analytic solutions for the long-range poten- 
tial |l6l Il7| , and the angular momentum insensitivity of 
a properly defined short-range K matrix K c EHIl. We 
show that as far as K c is concerned, the hyperfine in- 
teraction can be ig nored, and the frame transformation 
EH [13 E3l applies basically exactly. This con- 
clusion greatly simplifies the description of any atomic 
collision that involves hyperfine structures. In the case 
of a collision between any two alkali-metal atoms in their 
ground state, whether they are identical or not, it reduces 
a complex multichannel problem to two single channel 
problems. This property, along with the energy and 
angular-momentum insensitivity of K c |l3l Il5| , leads to 
an accurate and complete characterization of slow col- 
lisions between any two alkali-metal atoms, including 
shape resonances, Feshbach resonances, practically all 
partial waves of interest, and over an energy-range of 
hundreds of millikelvins, by four parameters for atoms 
with identical nuclei, and five parameters for different 
atoms or different isotopes of the same atom. To be more 
specific, the four parameters can be taken as the singlet s 
wave scattering length aos, the triplet s wave scattering 
length oot, the C§ coefficient for the long-range van der 
Waals potential —C^/r 6 , and the atomic hyperfine split- 
ting A£Jf F (The reduced mass /i, which is also needed, 
is not counted as a parameter since it is always fixed and 
well-known). For different atoms or different isotopes of 
the same atom, we need another hyperfine splitting for 
a total of five parameters. These results also prepare us 
for future analytic representations of multichannel cold 
collisions, when we restrict ourselves to a smaller range 
of energies. 



II. MAQDT 

An iV-channel, two-body problem can generally be de- 
scribed by a set of wave functions 
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Here $i are the channel functions describing all degrees 
of freedom other than the inter-particle distance r; and 
Fij (r) satisfies a set of close-coupling equations 

j j=i 

(2) 

where fi is the reduced mass; U is the relative angular 
momentum in channel i; E is the total energy; and Vy 
is the representation of inter-particle potential in the set 
of chosen channels (see, e.g., reference 9] for a diatomic 
system with hyperfme structures). 

Consider now a class of problems for which the poten- 
tial at large distances (r > ro) is of the form of 

V ij (r) = (E i -C ni /r ni )6 ij , (3) 

in the fragmentation channels that diagonalize the long- 
range interactions. Here ni > 2, and Ei is the threshold 
energy associated with a fragmentation channel i. As an 
example, for the scattering of two alkali-metal atoms in 
their ground state, the fragmentation channels in the ab- 
sence of any external magnetic field are characterized by 
the FF coupling of reference [jj; differences in thresh- 
old energies originate from atom hyperfme interaction; 
?i; = 6 corresponds to the van der Waals interaction; and 
ro, with an order of magnitude around 30 a.u., corre- 
sponds to the range of exchange interaction. 

Before enforcing the physical boundary condition 
(namely the condition that a wave function has to be 
finite everywhere) at infinity, Eqs. J2J have N linearly 
independent solutions that satisfy the boundary condi- 
tions at the origin. For r > ro, one set of these solutions 
can be written as 

N 
i=l 

Here ff and gf are the reference functions for the long- 
range potential, —C n Jr ni , in channel i, at energy ej = 
E — Ei. They are chosen such that they are independent 
of both the channel kinetic energy and the relative an- 
gular momentum ij at distances much smaller than the 
length scale (3 ni = (2/jC„ i /?i 2 ) 1 ^™ i_2 ^ associated with 
the long- rang e interaction (see Appendix ^ and refer- 
ences [lj.llq '). 

Equation defines the short-range K matrix K c . It 
has a dimension equal to the total number of channels, 
N, and encapsulates all the short-range physics. The K c 
matrix can either be obtained from numerical calcula- 
tions (see Appendix[BJ or be inferred from other physical 
quantities such as the singlet and the triplet scattering 
lengths, as discussed later in the article. 

At energies where all N channels are open, the solu- 
tions given by Eq. 10} already satisfy the physical bound- 
ary conditions at infinity. Using the asymptotic behav- 
iors of reference functions f c and g° at large r (see Ap- 
pendix ^ and reference jig), it is easy to show from 



Eq. (0} that the physical K matrix, defined by Eqs. (4) 
and (5) of reference is an Nx N matrix given in terms 
of K c by 

K(E) = -{Z% - Z° g K°)(Z h - Z^KT 1 ■ (5) 

Here Z c fg ,Z c gg ,Z c ff , and Z c gf are N x N diag- 
onal matrices with diagonal elements given by 
Zf g *\ei,li), Zgg *\ej, lj), Z *\ej, and Z g ^ (ei,k), 
respectively (see Appendix |A] and references [Til Il6^. 
Equation is of the same form as its single channel 
counterpart 0, 0] , except that the relevant quantities 
are now matrices, and K c is generally not diagonal. 

At energies where N Q of the channels are open (e^ > 0, 
for i G o), and N c = N — N Q of the channels are closed 
(e^ < 0, for i € c), the physical boundary conditions at 
infinity leads to N c conditions that reduce that number 
of linearly independent solutions to N Q 0, 0, [2l| . The 
asymptotic behavior of these N Q solutions gives the N x 
N Q physical K matrix 

K{E) = -{Z% - Z° gg Kl ff )(Z° n - ZfrKSf,)- 1 . (6) 

Here Zj g , Z gg , Zjj, and Z g * are N Q x N diagonal ma- 
trices with diagonal elements given by the corresponding 
Z c matrix element for all open channels; and we have 
defined the effective K c matrix for the open channels, 
K^f, to be 

Kl f S = K oo + K oc(x c - K c )- l Ko ■ (7) 

Here \ c i s an N c x N c diagonal matrix with elements 
X c ( n »)( £i) y (see Appendix 1X1 and references 0,0) for 
all closed channels. K£ , K£ c , K£ , K£ c , are submatrices 
of K c corresponding to open-open, open-closed, closed- 
open, and closed-closed channels, respectively. 

All on-the-energy-shell scattering properties can be de- 
rived from the physical K matrix. In particular, the 
physical S matrix is given by 

S(E) = [I + iK(E)][I-iK(E)]- 1 , (8) 

where I represents a unit matrix. From the S matrix, 
the scattering amplitudes, the differential cross sections, 
and other physical observables associated with scattering 
can be easily deduced Q. 

It is worth noting that Eq. © preserve the form of 
Eq. JSJ. Thus the effect of closed channels is simply to 
introduce an energy dependence, through x c , into the 
effective K c matrix, K^f, for the open channels. In par- 
ticular, the bare (unshifted) locations of Fcshbach reso- 
nances, if there are any, are determined by the solutions 
of 

det[ X c (E) - K c cc ] = . (9) 

They are locations of would-be bound states if the closed 
channels are not coupled to the open channels. The same 
equation also gives the bound spectrum of true bound 
states, at energies where all channels are closed. 
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This completes our summary of MAQDT. It is com- 
pletely rigorous with no approximations involved. The 
theory is easily incorporated into any numerical calcula- 
tions (see Appendix lB|l . The difference from the standard 
approach is that one matches the numerical wave func- 
tion to the solutions of the long-range potential to extract 
K°, instead of matching to the free-particle solutions to 
extract K directly. This procedure converges at a much 
smaller r = ro, the range of the exchange interaction, 
than methods that match to the free-particle solutions. 
Furthermore, since the propagation of the wave function 
from ro to infinity is done analytically, through the Z c 
matrix for open channels and x° function for closed chan- 
nels, there is no difficulty in treating shape resonances 
right at or very close to the threshold. This improved 
convergence and stability does not however fully illus- 
trate the power of MAQDT formulation and is not the 
focus of this article. Instead, we focus here on the simple 
parameterization of slow atomic collisions with hypcrfinc 
structures made possible by MAQDT. The result also 
lays the ground work for future analytic representations 
of cold collisions. 



III. SIMPLIFIED PARAMETERIZATION WITH 
FRAME TRANSFORMATION 

Equations JSJ-JTJl, and © already provide a param- 
eterization of slow-atom collisions and diatomic bound 
spectra in terms of the elements of the K c matrix. For 
alkali-metal atoms in their ground state, where the mul- 
tichannel nature arises from the hyperfine interaction, or 
a combination of hyperfine and Zeeman interactions for 
scattering in a magnetic field, this parameterization can 
be simplified much further by taking advantage of a frame 
transformation [1 fTo L 111 Im E ^ . 

At energies comparable to, or smaller than the atomic 
hyperfine and/or Zeeman splitting, one faces the di- 
chotomy that the hyperfine and/or Zeeman interaction, 
while weak compared to the typical atomic interaction 
energy, is sufficiently strong that the physical K matrix 
changes significantly over a hyperfine splitting. (This 
is reflected in the very existence of Feshbach resonances 
P, and states with binding energies comparable to or 
small than the hyperfine splitting.) As a result, the frame 
transformation does not apply directly to the physical K 
matrix itself, and is generally a bad approximation even 
for the K° matrix of reference |3 ■ It was this recogni- 
tion that first motivated the solutions for the long-range 
potentials |la. Il7|. 

This dichotomy is easily and automatically resolved 
with the introduction of the short range K matrix K c . 
The solution is simply to ignore the hyperfine and/or 
Zeeman interaction only at small distances and treat it 
exactly at large distances. For r < r , the atomic in- 
teraction is of the order of the typical electronic energy. 
Thus as far as K c , which converges at ro, is concerned, 
the hyperfine and/or Zeeman interaction can be safely 



ignored. In this approximation, the K c matrix in the 
fragmentation channels can be obtained from the K c ma- 
trix in the condensation channels, namely the channels 
that diagonalize the short-range interactions, by a frame 
transformation. 

For simplicity, we restrict ourselves here to the case 
of zero external magnetic field, although the theory can 
readily be generalized to include a magnetic field. The 
fragmentation channels are the FF coupled channels 
characterized by quantum numbers [t| : 

(a 1 L 1 S 1 J 1 I 1 F 1 ) A (a 2 L 2 S 2 J 2 I 2 F 2 ) B FlTM T , 

where F results from the coupling of F\ and F 2 ; I is 
the relative orbital angular momentum of the center-of- 
masses of the two atoms. T represents the total angular 
momentum, and Mt is its projection on a space-fixed 
axis 0- 

Provided that the off-diagonal second-order spin- 
orbital coupling |22j | can be ignored, a good approxima- 
tion for lighter alkali-metal atoms, or more generally, for 
any physical processes that are allowed by the exchange 
interaction, the condensation channels can be taken as 
the LS coupled channels characterized by quantum num- 
bers @: 

{a x L x Sxh)A{oL2L2S 2 h) B LlCSKITM T , 

where C = L + 1 is the total orbital angular momentum. 
S = Si + S2 is the total electron spin. I = T + 1 2 is the 
total nuclear spin. And K = C + S is the total angular 
momentum excluding nuclear spin. 

Ignoring hyperfine interactions, as argued earlier, the 
K c matrix in Fi^-coupled channels, labeled by index i or 
j, is related to the K c matrix in LiS-couplcd channels, 
labeled by index a or /3, by a frame transformation 



a/3 



K« LS >U j0 



(10) 



where K C ( LS ^ is the K c matrix computed in the LS cou- 
pling with the hyperfine interactions ignored. The most 
general form of frame transformation Ujp is given by 
Eq. (49) of reference 0|. For collision between any two 
atoms with zero orbital angular momentum, L± = L 2 = 
L = 0, including of course any two alkali-metal atoms in 
their ground states, the frame transformation simplifies 
to 

U lP {T) = 6 kh (-l) F >+ s e +I e[F li ,F 2i ,F i ,Sf 3 ,K /3 ,I f3 } 1 / 2 
P 1 T 1 I 1 2 & 



Fu F 2 i Fi 



for atoms with different nuclei. Here [a, b, . . .] = (2a + 
1) (26 + 1) For two atoms with identical nuclei, the 
same transformation needs to be multiplied by a normal- 
ization factor 

U {i}{P) = {^ + S{a2L 2 S 2l a 1 L 1 S 1 )[l-6(J 2t F 2l ,J ll F ll )}} 1 / 2 
xU ip . (12) 
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We emphasize that to the degree that the hyperfine in- 
teraction in a slow atomic collision can be approximated 
by atomic hyperfine interactions, as has always been as- 
sumed |2|, the frame transformation given by Eq. I|1U|) 
should be regarded as exact. If the hyperfine interaction 
inside ro, the range of the exchange interaction, cannot 
be ignored, the true molecular hyperfine interaction 
would have to be used. Inclusion of atomic hyperfine in- 
teractions inside ro is simply another approximation, and 
an unnecessary complication, that is of the same order of 
accuracy as ignoring it completely. In other words, any 
real improvement over the frame transformation has to 
require a better treatment of molecular hyperfine inter- 
actions |'2.'"{l| . A similar statement is also applicable to the 
Zeeman interaction. 

The applicability of the frame transformation greatly 
simplifies the description of any slow atomic collision 
with hyperfine structures. For alkali-metal atoms in 
their ground state, and ignoring off-diagonal second- 
order spin-orbital coupling |22j , it reduces a complex mul- 
tichannel problem to two single channel problems, one 
for the singlet S = 0, and one for the triplet, 5=1, 
with their respective single-channel K c 0, Q denoted 
by Kg{e,li) and K^(e,k), respectively. The K c matrix 
in the LS coupling, K C ^ LS \ is diagonal with diagonal ele- 
ments given by either Kg or |9j . Ignoring the energy 
and the angular momentum dependences of K$(e, U) and 
K^(e,k) [lJJ Uja], they become simply two parameters 
K% = #f (0,0) and K% = K£(0,0), which are related to 
the singlet and the triplet s wave scattering lengths by 

M 



■2b 



T(l-b) 



r(i + 6) 



FT c (0,0) + tan(7r&/2) 
if c (0,0)-tan(7rV2) 



(13) 



where b = l/(n — 2) with n = 6 for alkali- metal scat- 
tering in the ground state. With K C ^ LS \ and therefore 
K c , being parameterized by two parameters, a complete 
parameterization of alkali-metal scattering requires only 
two, or three, more parameters including Cq, which de- 
termined the length and energy scales for the long range 
interaction, and the atomic hyperfine splitting AE^ F , 
which characterizes the strength of atomic hyperfine in- 
teraction and also determines the channel energies. 

We note here that our formulation ignores the weak 
magnetic dipole-dipole interaction [22T l25j . It is im- 
portant only for processes, such as the dipolar relax- 
ation, that are not allowed by the exchange interaction. 
Such processes can be incorporated perturbatively after 
a MAQDT treatment [llj. We also note that for pro- 
cesses, such as the spin relaxation of Cs, for which the 
off-diagonal second-order spin-orbital coupling is impor- 
tant |22l |26]| . a different choice of condensation channels, 
similar to the J J-coupled channels of reference , would 
be required. The resulting description is similar concep- 
tually, but involves more parameters [2^, [2(| . 




E/k„ ( K ) 



FIG. 1: \Su - l\, where S a is an S matrix element, for 
the s wave elastic scattering of two 23 Na atoms in channel 
[{Fi = 1,F 2 = 1}F = 0,1 = 0,T = 0], as a function of 
F/fes, where fcs is the Boltzmann constant. The vertical 
lines identify the locations of thresholds for {Ft = l,Fa = 
2} and {Fi = 2, F2 = 2} channels. Solid line: results of a 
four-parameter MAQDT parameterization. Points: results of 
close-coupling calculations. 



IV. SAMPLE RESULTS FOR SODIUM-SODIUM 
SCATTERING 

As an example, Figures 11131 show the comparison be- 
tween close-coupling calculations and a four-parameter 
MAQDT parameterization for slow atomic collision be- 
tween a pair of 23 Na atoms in the absence of exteranl 
magnetic field. The points are the close-coupling results 
using the potentials of references [13, 13 ■ The curves 
represent the results of a four-parameter parameteriza- 
tion with aos = 19.69 a.u., oqt = 64.57 a.u., Cq = 1556 
a.u. [2cJ, and SE^ F = 1772 MHz, where a os and a 0T 
are computed from the singlet and the triplet poten- 
tials of references [23, H5|- Figure 1 shows the S ma- 
trix element for the s wave elastic scattering in chan- 
nel [{Fi = 1,F 2 = 1}F = 0,1 = 0,T = 0]. The fea- 
ture around 130 mK is a Feshbach resonance in channel 
[{Fi = 2, F 2 = 2}F = 0, I = 0,T = 0]. For this particular 
case, K C ( LS ) is a 2 x 2 matrix 



R c(LS) 



K% 
KZ 



(14) 



with channel ordering shown in Table [I] Kg and K^ are 
related to the singlet and the triplet scattering lengths 
by Eq. (|13f) . The frame transformation is given by [c.f. 
Eqs. (HU and JX2J] 



U(T = 0) = 



1 



V3 



2V2 \ VZ — v/3 



which leads to 



K c = 



?>K% + 5KZ 



VS{K%-K^ 



8 I VT5 {K% - Kf) 5K C S + 3K° 



(15) 



(16) 
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TABLE I: Channel structure for s wave scattering between 
two identical atoms with L 1 = L2 = 0, Si = S2 = 1/2, and 
Ii = 1% = 3/2, in the absence of external magnetic field. 
Examples include 7 Li, 23 Na, 39 K and 87 Rb. 



T 


LS coupling (S,I) 


FF coupling {Fi,F 2 }F 





S=0, 1=0 


{1,1}0 




S=l, 1=1 


{2,2}0 


1 


S=l, 1=1 


{1,2}1 


2 


S=0, 1=2 


{1,1}2 




S=l, 1=1 


{1,2}2 




S=l, 1=3 


{2,2}2 


3 


S=l, 1=3 


{1,2}3 


4 


S=l, 1=3 


{2,2}4 




E/k„ ( K ) 



FIG. 2: The same as Figure 1 excepts it is for d wave channel 
[{Fi = 1,F 2 = 1}F = 2,1 = 2,T = 2]. 




E/k B ( K ) 



FIG. 3: The 5* matrix element, |iSy| , for the s wave inelastic 
scattering of two 23 Na atoms between channels [{Fi = 1, F2 = 
1}F = 2, 1 = 0, T = 2] and [{Fi = 1, F2 = 2}F = 2, 1 = 0, T = 
2]. 



because the [{Fi = 1, F 2 = 1}F = 0, / = 0, T = 0] chan- 
nel is not coupled to {F± — 1, F 2 = 2} channels. 

The agreements between the MAQDT parameteriza- 
tion and close-coupling calculations are excellent, exact 
for all practical purposes, in all cases. Conceptually, 
these results illustrate that through a proper MAQDT 
formulation, atomic collision over a wide range of en- 
ergies (300 mK compared to the Doppler cooling limit 
of about 0.2 mK for 23 Na), with complex structures in- 
cluding Fcshbach and shape resonances, and for different 
partial waves, can all be described by parameters that 
we often associate with the s wave scattering at zero en- 
ergy only, namely the singlet and the triplet scattering 
lengths. 



From the K c matrix, the S matrix is obtained from the 
MAQDT equations Note how Eq. (JTSJ shows 

explicitly that the off-diagonal element of K c , which de- 
termines the rate of inelastic collision due to exchange 
interaction, goes to zero for K s = K^, namely when 
aos = &0T- 

The results presented in Figs. 2 and 3 are obtained 
in similar fashion. Figure 2 shows the S matrix ele- 
ment for the d wave elastic scattering in channel [{Fi = 
1,F 2 = 1}F = 2,1 = 2,T = 2]. It illustrates how the 
same parameters that we use to describe the s wave 
scattering also describe the d wave scattering, due to 
the fact that Kg and K^ are insensitive to I [l3l [l5| . 
Here the sharp features around the thresholds are d wave 
shape resonances. Figure 3 shows the S matrix ele- 
ment for the s wave inelastic scattering between chan- 
nel [{Fi = 1,F 2 = 1}F = 2,1 = 0,T = 2] and channel 
[{Fi = 1,F 2 = 2}F = 2,1 = 0,T = 2]. The kinks (dis- 
continuities in the derivative), in both Fig. Eland Fig.Q] 
at the {Fi = 2, F 2 = 2} threshold, are general features 
associated with the opening of an s wave channel. There 
is no kink at the {F 1 = 1,F 2 = 2} threshold in Figure □ 



V. CONCLUSION 

In conclusion, a multichannel, angular-momentum- 
insensitive, quantum-defect theory (MAQDT) for slow 
atomic collisions has been presented. We believe it to 
be the optimal formulation for purposes including ex- 
act numerical calculation, parameterization, and ana- 
lytic representation. We have shown that by dealing 
with the short-range K matrix K € , the frame transfor- 
mation becomes basically exact, which greatly simplifies 
the description of any slow atomic collision with hyper- 
fine structures. As an example, we have shown that 
even a simplest parameterization with four parameters, 
in which the energy and the I dependence of K s and 
Fy are completely ignored, reproduces the close-coupling 
calculations for 23 Na atoms over a wide range of energies 
basically exactly. The effect of an external magnetic field, 
which is not considered in this article, is easily incorpo- 
rated as it simply requires another frame transformation 

The concepts and the main constructs of the theory 
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can be generalized to other scattering processes including 
ion-atom scattering and atom-atom scattering in excited 
states. The key difference will be in the long-range in- 
teraction [c.f. Eq. J2J)]- In addition to possibly different 
long-range exponent (such as rij = 4 for ion-atom scat- 
tering), there may also be long-range off-diagonal cou- 
pling that will have to be treated differently. 

Finally, we expect that if we restrict ourselves to a 
smaller range of energies, of the order of (h 2 /2fi)(l/ fo) 2 
(about 1 mK for 23 Na), a number of analytic results, sim- 
ilar to the single-channel results of references 1 and @, 
can be derived even for the complex multichannel prob- 
lem of alkali-metal collisions. These results may, in par- 
ticular, lead to a more general and more rigorous param- 
eterization of magnetic Feshbach resonances (see, e.g., 
references 12, 30] for some recent works in this area). 
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which can be written in a dimensionless form as 



d 2 _ 1(1 + 1) 
dr 2 r? 



u ts i{r s ) = 



(A2) 



where r s — r/(3 n is a scaled radius, j3 n = 
(2/j,C n /h 2 ) 1 ^ n ~ 2S> is the length scale associated with the 
—C n /r n interaction, and 

(A3) 



{V/2n){l/p n y 

is a scaled energy. 

The f c and g c pair are chosen such that they have not 
only energy-independent, but also angular-momentum- 
independent behaviors in the region of r <C (3 n (namely 
r s < 1): 



fli(r s ) 



r»<l 



r»<l 



(2/7r) 1/2 r™/ 4 cos(y-7r/4) , (A4) 
-(2/ 7 r) 1 / 2 r^ 4 sin( 2/ -V4) , (A5) 



APPENDIX A: DEFINITIONS OF MAQDT 
FUNCTIONS 



for all energies [H 0. Here y = [2/(n - 2)]r s ( " 2)12 . 
They are normalized such that 



The reference functions f c and g c for a —C n /r n [n > 2) 
potential are a pair of linearly independent solutions of 
the radial Schrodinger equation 



~2ndr 2 + 



h 2 l(l + 1) 
2/xr 2 



u d (r)=0, (Al) 



W(f c , g c ) = r A -f - - 2/tt . (A6) 

dr s dr a 

For e = 0, the f c and g c pair for arbitrary n can be 
found in reference |15| . For e ^ 0, the f c and g c pair for 
n = 6 can be found in referenc e l3ll . The are related to 
the f° and g° pair of reference [16| by 



l c \ J_ f cos(7r^ /2) sin(™ /2) 1 \ I f° \ 
g c V2 I - am(nv /2) cos(™ /2) ) \ -1 j \ g° j ' 



(A7) 



where vq = (21 + l)/4 for n = 6. 

The Z c ( n \e s , I) matrix is defined by the large-r asymptotic behaviors of f c and g c for e > 



J ff 



(n) 
f 

c(n) 



In 
~2 



7 c(n), 



[z^( £s ,0sin (k s r s - - Z«;He s ,l)cos ( 





ZttV 




~ v. 




ZttV 







(A8) 
(A9) 



where k s = k/3 n with k = (2^ie/h 2 ) 1 ^ 2 . This define a 2 x 2 
Z c{n \e si l) matrix 



The X^ n ^( c a) function is defined through the large-r 
asymptotic behaviors of f c and g c for e < 0. 



Z ff Z f ? 

At 



It is normalized such that 



9<J 



(A10) 



det 



7c(n) 



7 c(n) yc{n) 



ryc{n) ryc(n) 

Z 9S Z /9 



(All) 



7 



9eA r s) 



r c{n) 



c(n) . 



(27TK, 



v-1/2 



(«), 



(A12) 
(A13) 



r 



where n s = n[} n with k = (2i±\Ei\/h 2 ) 1 / 2 . This defines a 
2x2 W c W(e s ,0 matrix, 



from which the Xi^ n \ e a) function is defined by 

c(rl) / \ Ti^cfn) /jrrc(n) 

Xl ( >(e s ) = W f l J /Wgi • 
The W c ^ matrix is normalized such that 



w c(n) w c(n) 
w c(n) w c(n) 



(A14) 



(A15) 



dct 



w c(n) c(„) _ w o K n >w c^, = j _ (Alg) 



c (™)w c ( n ) 



The Z c (")(e s , and W c ( n )(e a ,f) matrices, for e > and 
e < 0, respectively, describe the propagation of a wave 
function in a — C n /r n potential from small to large dis- 
tances, or vice versa. They are universal functions of the 
scaled energy e s with their functional forms determined 
only by the exponent n of the long-range potential and 
the I quantum number. The C n coefficient and the re- 
duced mass play a role only in determining the length 
and energy scales. 

The Z c{ ^ matrix for n = 6 is given in reference (l4| . 

The Xt \ e s) function for n = 6 is given in reference 
[l3| . They are derived from Eq. (|A7|I and the asymptotic 
behaviors of the f° and g° pair given in reference . 



APPENDIX B: K c FROM NUMERICAL 
SOLUTIONS 

Let F(r) be the matrix, with elements Fij(r), repre- 
senting any N linearly independent solutions of the close- 



coupling equation, and F'(r) be its corresponding deriva- 
tive [Each column of F(r) corresponds to one solution 
through Eq. (J)] . For r > r , F can always be written as 



F(r) =f c {r)A-g%r)B 



(Bl) 



where f°(r) and g c {r) are N x N diagonal matrices with 
diagonal elements given by ff{r) and gf(r), respectively. 
The matrices A and B can be obtain, e.g., from knowing 
F(r) and F'(r) at one particular r > r$. Specifically 



.4 
B 



(7r(3 n /2)[g c '(r)F(r) 
(7r/3„/2)[/ c '(r)f(r) 



g c (r)F'(r)} . 
- f C (r)F'(r)} 



(B2) 
(B3) 



Comparing Eq. (|B1|I with Eq. J2J gives 



K c = [/°'(r)F(r)-/ c (r)F'(r)]|5 c '(r)F(r)- ff c (r)F / (r)]- 1 . 

In an actual numerical calculation, which can be imple- 
mented using a number of different methods [33, the 
right-hand-side (RHS) of this equation is evaluated at 
progressively greater r until K c converges to a constant 
matrix to a desired accuracy. This procedure also pro- 
vides a numerical definition of , namely it is the radius 
at which the RHS of Ea. (|B4|) becomes a r-independent 
constant matrix. 
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